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Abstract 

A theoretical study of laser and plasma based electron acceleration is presented. An effective model has been used, in 
which the presence of an underdense plasma has been taken account via its index of refraction n^. In the confines of 
this model, the basic phenomena can be studied by numerically solving the classical relativistic equations of motion. 
The key idea of this paper is the application of chirped, bichromatic laser fields. We investigated the advantages and 
disadvantages of mixing the second harmonic to the original A = 800 nm wavelength pulse. We performed calculations 
both for plane wave and Gaussian pulses. 

Keywords: Underdense plasma. Electron acceleration. Classical electrodynamics. Relativistic equation of motion. 
Ultrashort laser pulses 


1. Introduction 

The history of laser-plasma based electron acceleration 
began in the late ’70-s. Tajima and Dawson predicted 
that the plasma wakes, generated by the ponderomotive 
force of short, intense laser pulses, are capable to accelerate 
bunches of electrons effectively [1]. In this scheme, there 
are many methods for generating high amplitude plasma 
wakes, for a summary, see Ref. [2] . The key point in every 
method is the resonant excitation of the plasma. In the 
middle ’80-s, the invention of the Chirped Pulse Amplifi¬ 
cation (CPA) made it possible to generate short, intense 
laser pulses without damaging the medium [3]. This new 
technology is an important milestone both in the history 
of laser physics and compact, plasma based particle accel¬ 
erators. Nowadays, thanks to the advanced technological 
developments, the vision of Tajima and Dawson is getting 
reality. 

There is a serious demand for that, since the conven¬ 
tional storage ring technology has reached its limits: the 
amplitude of the accelerating gradient must not be larger 
than 50 MV/m. The violation of this criterion would lead 
to electric discharges that would critically damage the sys¬ 
tem. The present state-of-the-art technology, namely, the 
CERN-LHC, is based on 8.3 T strong superconducting 
magnets. This way, the theoretical maximum of center- 
of-mass (CM) energy is 14TeV. According to the most re¬ 
cent news, the LHC is operating at 13TeV at the moment. 
Using the standard technology, the maximal CM energy 
can be improved in two independent ways. One either 
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has to build a larger ring (VLHC), with a circumference 
of 80 or 100 kms, or develop stronger (e.g. 21T) super¬ 
conducting magnets. Both solutions would be extremely 
expensive, and the latter is also very uncertain, there is no 
guarantee for it to succeed. Due to theses difficulties, new 
technologies are needed. The most popular of them is the 
concept of laser-plasma based particle accelerators, that 
has been mentioned above. Nowadays, there are promis¬ 
ing experimental results for building compact particle ac¬ 
celerators: electrons have been accelerated up to multiples 
of GeVs within a few cm long plasma cell [4, 5]. CERN 
is also open for new technologies: the construction of the 
CERN AWAKE experiment (Proton Driven Plasma Wake¬ 
field Acceleration) has already begun. The details of the 
proposed scheme can be found in [6, 7]. It is important to 
mention that this scheme has been designed for electron 
acceleration only. For a long time, it was thought that 
acceleration of positrons by plasma wakes is impossible. 
Recently, it has been shown that by applying a “doughnut 
shaped” driver pulse, laser-plasma based positron acceler¬ 
ation can be realized as well [8]. 

This study is the sequel of our recent work [9]. In that 
paper, an effective theory for describing electron accel¬ 
eration in underdense plasmas has been presented. We 
showed that a single electron can be effectively accelerated 
both by monochromatic planewave pulses and Gaussian 
laser pulses, up to 270 MeV, which agrees quite well with 
other theoretical and experimental results. In the present 
paper we investigate the advantages of applying a bichro¬ 
matic driver pulse, namely, adding the second harmonic to 
the original laser pulse. Ehlotzky’s work [10], which sum¬ 
marizes various relevant atomic phenomena in bichromatic 
laser fields gave us a good motivation and serves as a start¬ 
ing point of our present work. To our knowledge, there 
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is no such approach in the literature. However, there is a 
nice proposal for producing narrow-energy-spread electron 
bunches from laser wakefield acceleration, using bichro- 
matic laser pulses [11]. 

In Section 2, the fundamental theoretical basics of laser- 
plasma based electron acceleration are overviewed, within 
the confines of our effective theory. The results are pre¬ 
sented in Section 3. Finally, we summarize our work. 


2. Theory 


During the laser-electron interaction, the Lorentz-force 
drives the motion of the electron: 


F = e (E -I- V X B) 


( 1 ) 


with e the electron charge, E the electric field, B the 
magnetic field, v the velocity of the electron and F the 
Lorentz-force. At sufficiently high intensities, the electron 
becomes relativistic. Therefore, one has to solve the rela¬ 
tivistic Newton-Lorentz equation: 


dp 

dt 


= e 




dy 1 
dt rrieC^ 


(2a) 

(2b) 


It is known that the electromagnetic field has to satisfy 
the electromagnetic wave equation. This condition yields 
the most general form for the electric and magnetic field: 


E(t,r) =££'o/[w0 (t,r)], 

(3) 

B(t,r) = -n X E(t) 

(4) 


with £ the polarization vector, Eq, the amplitude of the 
electric field, oj the angular frequency and n the unit vec¬ 
tor of the propagation of the electromagnetic field, respec¬ 
tively. / may be any arbitrary, smooth function. For 
a better transparency, we introduced the following nota¬ 
tion, since the electromagnetic field depends only on the 
planewave-argument 

0(t,r) := f — n • (5) 

c 


which depends on the laser and plasma frequencies lol and 
u}p, respectively, in the following way: 


— 



(7) 


with 


Some 


( 8 ) 


with rig the electron density in the plasma, Eq the permit¬ 
tivity of vacuum and rrie the electron mass. 

This approximation yields an effective theory which is 
useful for performing basic studies of the laser and plasma 
based electron acceleration via numerically solving the rel¬ 
ativistic equations of motion (2), treating E and B as a 
function of Q{t,r,nm)- The = 1 case describes the 
pure laser based electron acceleration. It is important to 
mention that a single electron cannot gain a net energy 
from a plane wave pulse in vacuum since during one pulse 
period, the electron gains and loses the same amount of 
energy during oscillating in the laser field. However, this 
symmetry can be broken by applying a (linear) chirp to 
the laser frequency: 


uj{t) = Wo + o-t, 


(9) 


(T being the chirp parameter and wq the initial frequency 
of the laser. 

It is convenient and useful to rescale the parameters and 
introduce the following dimensionless variables: 

cEq 

oo = -, 

nieUjQC 

Wot —> t, —r —> r, -p, (10) 

c rrieC 

wq © — y ©, —^ ^ (T- 

w^ 

The equations of motion, expressed in terms of these new 
variables, take the following form: 


dp 



dy 

dt 


do 

7 


E • p. 


(lla) 

(llb) 


If we also want to take into account the presence of 
a medium with an index of refraction < 1—for in¬ 
stance, an underdense plasma, which is the field of our 
present investigation—, we need to generalize the defini¬ 
tion of 0(t,r) in the following way: 

r 

e(t,r,n^) :=t-nmn- (6) 

c 

We interpret the generalized definition of 0 such that it de¬ 
scribes the electron propagation in an underdense plasma 
[12]. All the background effects are incorporated into Um, 


In general, a bichromatic electromagnetic field can be 
expressed in the compact form of 

A . , 

E — El H-Eg, (12a) 

B = Bi -b -Bg (12b) 

9 

with El and Bi being the electric and magnetic fields of 
the main harmonic and 0 < A < 1 the relative amplitude 
of the harmonics, q denotes the index of the (higher) har¬ 
monic with gwo initial frequency. The q~^ factor is a direct 
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consequence of the definition of the intensity parameter oq 
( see Eq. (10)): 

a.(«) = (13) 

For first, we define the harmonic of a bichromatic 
plane wave pulse with a sine-square shaped temporal en¬ 
velope: 

f (<S) ^ 

/ (0) = < sin (^qQ + q'^aqQ^ + (fq) if0S[O,T] (^4) 
[ 0 otherwise 

with T the pulse duration, aq the dimensionless chirp pa¬ 
rameter and ipq with q = 1 the carrier-envelope phase and 
ipq with q > 1 the relative phase of the harmonics. The 
electric field is polarized in the x direction and propagates 
in the y direction, that is, e = e^,, n = ej,. 

Gaussian pulse shapes provide a more realistic descrip¬ 
tion of laser beams. The mathematical expressions for an 
a;-polarized Gaussian beam that propagates in the 0 di¬ 
rection can be derived from the paraxial approximation 
[13, 14]. The explicit expressions containing first order 
corrections can be found in many works, see e.g. [15, 16]. 
Here we only present the generalizations of these formulae 
for the ( 7 *^ harmonic in terms of dimensionless variables: 


l^q,x = T7777\ 


Wq{z) 


COS 


W^{z) 


exp 


02 


(woT)^ 


2Rq{z) 


- ^q{z) + qQ + q^aQ"^ + (pg 


Eq^y = 0, 


(15a) 


(15b) 
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Rq{zr^ 
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_ lEo 
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sin 
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(oJoT) 

(15c) 


and the magnetic field is given by 


Rq{z) = (a;o/c) 2 : [l-I-(z,^/j/z)2] the radius of curvature, 
^q{z) = ta.n~^{z/Zq^B.) the Guoy phase, Zq^a the Rayleigh 
length, kq = g(l -|- qaqQ) the dimensionless wavenumber 
of the harmonic, Aq the initial wavelength of the main 
harmonic and T the pulse duration. 

The construction of a bichromatic Gaussian pulse im¬ 
plies some problems. Theoretically, there is a freedom 
by defining the parameters of the individual harmonics: 
one has to choose either the beam waists or the Rayleigh 
lengths to be equal. Since the beam waist and the Rayleigh 
length depend on each other, equal Rayleigh lengths would 
result in different beam waists. From the expressions given 
above if follows that the main and the higher harmonics 
have the same beam waist, this is the physically relevant 
choice. However, the equality of the beam waists implies 
the difference of the Rayleigh lengths. It is also clear from 
equations (15), (16) and (12) that the precise mathemat¬ 
ical expression of a bichromatic Gaussian pulse is rather 
complicated. 


3. Results 

In this section we discuss the advantages and disadvan¬ 
tages of applying bichromatic (g = 2) laser pulses for laser 
and plasma based electron acceleration. As discussed in 
our earlier work [9] , within the confines of the present ap¬ 
proach the Um < 1 case agrees quite well with the rim = 1 
case, since at relevant plasma densities the refraction index 
of an underdense plasma differs negligibly from unity. 

We investigated the effect of the presence of the second 
harmonic. We found that by properly chosen parameters 
the energy gain of the electron can be enhanced by 4 % or 
more. Graphically said, this is due to the increment of the 
accelerating field at the sharp rising and falling edges at 
the last oscillations of the electric field (see Fig. 1). 

In our calculations we chose the relative phase of the two 
harmonics to be zero. We compared the energy gain pro¬ 
vided by a monochromatic and a bichromatic pulse such 
that the two pulses have the same intensity. This means 
that if the intensity parameter of a monochromatic pulse 
is ao, then the intensity parameter of the corresponding 
bichromatic pulse is bao with 


Bq,x = 0, 

Eq^y — Eqx-, 

B„ j, = — ^^——E„ x~ 
Rq{z) 


(16a) 

(16b) 
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with Wo = {ojQ/c) ^/XoZqJi/n the beam waist, Wq{z) = 
(wo/c) \l + {z/Zq^nY]^ the beam radius at distance z. 


\ 1/2 

|Ei|^d0 


/•©/ 

A 

2 

/ 

El -k —E 2 

d0 

\-JQi 

2 

/ 


with 0i and 0/ being the plane wave arguments at times 
ti and tf as the interaction between the electron and the 
pulse starts and finishes. The subscripts i and / de¬ 
note “initial” and “final”, respectively. This normaliza¬ 
tion guarantees that the monochromatic and bichromatic 
pulses have the same intensity and validates the compari¬ 
son of the energy gains from the two different pulses. 
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The additional energy gain via the presence of the sec¬ 
ond harmonic had been investigated in the following way. 
We took a monochromatic pulse with a fixed pulse dura¬ 
tion at a given intensity and determined the optimal laser 
parameters and initial conditions that provide the most 
energy gain. After that we sought the optimal intensity 
ratio between the main and second harmonic while all the 
other parameters stayed fixed. At the determination of 
the optimal parameters we used the popular Nelder-Mead 
(also called “downhill simplex”) method [17]. 

[kV/nm] 



Figure 1: The x component of the electric field of a chirped, 
monochromatic (solid line) and bichromatic (dashed line) plane 
wave pulse. A = 800 nm, T = 35 fs, ao = 0.22, A = 0.48, 
CTi = —5.510 ■ 10“^, (T 2 = —1.311 ■ 10“^, (/?! = 0. Note the en¬ 
hancement of the accelerating field for the bichromatic case. 


At first we present the results for plane wave pulses, 
then for Gaussian pulses. We emphasize only the addi¬ 
tional energy gain caused by the presence of the second 
harmonic. That is, we normalize the energy gains (AA) 
to the optimal energy gain achieved by the corresponding 
a monochromatic pulse (AEq (A)) at the same intensity. 
The relative energy gain is defined by AE/AEq{A). We 
also investigate the behavior of the energy gain as a func¬ 
tion of the parameters of the two harmonics at fixed inten¬ 
sity ratio, namely, the chirp parameters and the carrier- 
envelope phase. In these cases we also scale the energy 
gain to unity, the normalization factor is the maximal en¬ 
ergy gain (denoted by AAmax) in a given parameter range. 
The scaled energy gain is also referenced as relative en¬ 
ergy gain as is defined by AE/AEynax- To avoid ambi¬ 
guity, the labels on the figures have been explicitly de¬ 
noted. For a better transparency we note that the values 
of ao = 0.22, 7, 12 and 22 correspond to the intensities of 
I = 10^^, 10^°, 3 • 10^° and 10^^ W/cm^, respectively. 

3.1. Plane wave pulses 

We found that the net energy gain depends very weakly 
on the chirp parameter of the second harmonic. The en¬ 
ergy gain is dominated by the chirp parameter of the main 
harmonic, this is indicated by the vertical stripes on Fig. 2. 
Practically this means that (J 2 can be chosen to be zero. 
As a consequence, both the experimental realization and 
numerical calculations are easier. 



Figure 2: The relative energy gain as a function of the dimensionless 
chirp parameters. Note that the <72-dependence is negligible. A = 
800 nm, T = 75 fs, ao = 22, A = 0.7, cpi = 0. 

AE 

-—[ 1 ] 

^^max 



Figure 3: The relative energy gain as a function of the carrier- 
envelope phase. A = SOOnm, T = 75fs, A = 0.7, ai = —5.518-10“^, 
(T 2 = —1.732 ■ 10“^, ao = 22 (solid line), ao = 7 (dashed line). Note 
the sensitivity to ifi\ 



Figure 4: The relative energy gain as a function of the relative am¬ 
plitude of the harmonics. The presence of the second harmonic 
enhanced the energy gain by about 4%. A = 800 nm, T = 75 fs, 
ao = 22, (T 1 = -5.530 ■ IQ-^, 02 = -1.732 ■ lO'^, ipi = 0. 
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Figure 5: The relative energy gain as a function of the carrier- 
envelope phase. A = 800nm, T = 5fs, ci = 8.494 ■ 10“^, (T 2 = 0, 
A = 0.24. ao = 7, Wo = 22.447r (solid line), ao = 12, Wq = 207r 
(dashed line). 
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Figure 6: The relative energy gain as a function of the chirp param¬ 
eters. AE depends very weakly on (T 2 . A = 800nm, T = 5fs, ao = 2, 
Wo = 207r, A = 0.1, = 2.058. 


It is not surprising that the carrier-envelope phase plays 
a very important role. However, it is not obvious that the 
energy gain depends very sensitively on there are two 
“worst-case” values at which the net energy gain is zero. 
There are also two optima, near 7 r /4 and 37 r /4 (see Fig. 3). 

Fig. 4 shows the key point of our study: at properly 
chosen laser parameters the energy gain of a single elec¬ 
tron can be enhanced with approximately 4 % by mixing 
the second harmonic to the main harmonic with a suit¬ 
able intensity ratio. In the next subsection it can also be 
seen that if the intensity ratio has not been set correctly, 
then the application of the second harmonic results in a 
net energy loss compared to the monochromatic case (see 
Fig. 7). 


AE 

-[ 1 ] 

AEo(A) 



3.2. Gaussian pulses 

In general, a Gaussian pulse can transfer the most en¬ 
ergy to a single electron if the electron initially moves on- 
axis an co-propagates with the beam. Namely, the initial 
momentum of the electron has the form of po = Po^z- 
This choice guarantees that the interaction length will be 
as high as possible. It is also important to initially place 
the electron far enough to not to feel the electric field of 
the pulse. The initial position of the electron has the form 
of xo = xqBz. The dimensionless initial position is spec¬ 
ified by Xo = ttuiq'T that corresponds to a ttcT distance 
from the center of the pulse. After these considerations it 
should be clear that the larger the value of po the larger 
the energy gain. This has been also confirmed by our nu¬ 
merical calculations, therefore we set po = 2 every time. 
This initial condition corresponds to 7 = \/5. 

For dichromatic Gaussian pulses we found that the ap¬ 
plication of the second harmonic is only reasonable at short 
pulse durations, namely in the 5fs — 15fs range. We anal¬ 
ysed the energy of the electron as a function of time dur¬ 
ing the interaction with the beam and found that at usual 
pulse durations, e.g. 35fs — 40 fs or above, the second har- 


Figure 7: The relative energy gain as a function of the relative am¬ 
plitude of the harmonics. The second harmonic enhanced the en¬ 
ergy gain by about 20 — 30 %, compared to the monochromatic case. 
A = 800nm, T = 5fs, cri = 8.494 ■ IQ-^, (T 2 = 0. oq = 7, Wo = 20n, 
ifii = 0 (solid line), ao = 12, Wo = 22.447r, tpi = 0 (dashed line). 


monic only causes a small, oscillatory perturbation in the 
energy-time function. 

However, at short pulse lengths, we found considerable 
additional gains due to the presence of the second har¬ 
monic. As expected, the energy gain depends very sensi¬ 
tively on the carrier-envelope phase. For some values of 
ipi, the net energy gain can be even negative (see Fig. 5). 

The energy gain depends very weakly on the chirp pa¬ 
rameter of the second harmonic, as presented on Fig. 6 , 
the dominant parameter—along with the initial momen¬ 
tum, carrier-envelope phase and beam waist—is the chirp 
parameter of the main harmonic. 

As mentioned above, we found that the application of 
the second harmonic may result in a considerable ad¬ 
ditional energy gain. The enhancement may be 20%, 
greater than conjectured, but it may reach even 30 % (see 
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Fig. 7). This promising result suggests that bichromatic 
laser pulses could be efficiently used for laser and laser- 
plasma based electron acceleration. 

4. Summary 

An effective theory for describing laser-plasma based 
electron acceleration has been presented. Earlier we 
showed that the background effects of the plasma can be 
incorporated into its refraction index. This way, the laser- 
driven plasma based acceleration can be well approximated 
with the pure laser-based acceleration, and the basic phe¬ 
nomena can be studied in a numerically and theoretically 
simple manner. In the present paper we investigated the 
acceleration mechanisms driven by bichromatic laser fields. 
We found that by properly chosen parameters the energy 
gain can be enhanced by 4% for plane wave pulses and 
even 30 % for Gaussian pulses, compared to the monochro¬ 
matic case. These are promising results that confirm that 
it is useful to apply bichromatic driver pulses for laser and 
laser-plasma based acceleration. It would be interesting 
to perform the same calculations with such Gaussian pulse 
shapes that are exact solutions of Maxwell’s equations and 
compare them with our most recent results. P. Varga and 
P. Torok derived such solutions and found that at wide 
focusing—that is, if the beam diameter is greater than ten 
times the laser wavelength—the paraxial approximation 
provides satisfactory results, however, at tight focusing, 
the differences are significant [18]. 
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